function diff = int_rep1_fun(m,param)

    % Alternative
    % n=10^5;
    % mm=linspace(param.m_h,m,n);
    % 
    % int=(param.zeta+delta_m(mm,param))./mm.^(1+param.gamma_R_1);
    % 
    % diff = sum(int./n).*(m-param.m_h);

    diff = nan(size(m));
    f = @(x) ( (param.zeta+delta_m(x,param))./x.^(1+param.gamma_R_1));

    for i=1:length(m)
        diff(i)=integral(f,param.m_h,m(i));
    end

end






